Snake States in Graphene p-n Junctions 



J. R. Williams^ and C. M. Marcus^ 

^School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA 
^Department of Physics, Harvard University, Cambridge, MA 02138, USA 
(Dated: December 8, 2010) 

We investigate transport in locally-gated graphene devices, where carriers are injected and col- 
lected along, rather than across, the gate edge. Tuning densities into the p-n regime significantly 
reduces resistance along the p-n interface, while resistance across the interface increases. This pro- 
vides an experimental signature of snake states, which zig-zag along the p-n interface and remain 
stable as applied perpendicular magnetic field approaches zero. Snake states appear as a peak in 
transverse resistance measured along the p-n interface. The generic role of snake states disordered 
graphene is also discussed. 



The low-energy spectrum of graphene, a two- 
dimensional hexagonal lattice of carbon atoms, yields two 
gapless modes of massless charge carriers that can be 
doped either n-type (electron-like) or p-type (hole-like) 
by contacts, electrostatic gates, or charged impurities 
on the graphene surface or within the supporting sub- 
strate. The unique band structure of graphene has sur- 
prising consequences in transport, including half-integer 
quantized conductance of 4(n + 1/2) e^//i in the quan- 
tum Hall regime [IHl], and finite minimum conductivity 
c"min ^ /h [11 [5]. Recently, the ability to locally con- 
trol density and carrier type by electrostatic gates has re- 
sulted in device configurations with adjacent n-type and 
p-type regions, separated by an electrically tunable p-n 
junction (PNJ) [6^^. The absence of a gap allows ballis- 
tic carriers that approach the PNJ normal to the interface 
to pass through without reflection, while carriers incident 
at oblique angles are reflected [9 . 

Besides intentionally gated p-n devices [SHH], p-n inter- 
faces play an important role conduction near the charge- 
neutrality point (CNP) in disordered graphene [10til2]. 
Near the CNP, disordered graphene consists of puddles 
of p and n regions whose boundaries — contours of zero 
density — form a percolating network of PNJs. Transport 
modes moving along p-n interfaces contribute to conduc- 
tion in the percolation region and alter the expected value 
of conductance, (Jmim at the CNP. 

In this Letter, we report low-temperature transport 
in a locally-gated graphene sample with one edge of 
the gate running between electrical contacts, allowing 
transport measurements along the PNJ [see Fig. 1(a)]. 
Transport is investigated along and away from the gate 
edge as a function of top gate voltage, Vt, back gate 
voltage, Fb, and perpendicular magnetic field, B. At 
zero and low magnetic fields, creating a PNJ along the 
gate edge decreases the longitudinal resistance along the 
edge and simultaneously increases resistance across the 
PNJ. These observations indicate the presence of cur- 
rent carrying modes confined to the p-n interface. As 
the perpendicular field is increased toward the quantum 
Hall regime, PNJ modes evolve into counter-propagating 
quantum Hall edge states moving along the p-n interface. 



with snake states altering transport at the transition be- 
tween hall plateaus. At low magnetic fields, increasing 
the electric field across the junction via gate voltages re- 
duces the contribution to conduction of the PNJ inter- 
face state, presumably reflecting a destabilization of edge 
modes with increased potential gradients. 

The basic mechanism leading to conduction along the 
PNJ is easily visualized at small but nonzero perpendicu- 
lar magnetic field, where a change in sign of charge carri- 




FIG. 1: (a) Schematic of the device. Contacts A and B are par- 
tially under and partially outside the top gate. Under certain values 
of Vt and V^, a PNJ forms (black dashed line), connecting con- 
tacts A and B. A close-up of the PNJ shows the modulation of the 
density across the junction that allows for an additional conduction 
channel to appear between contacts A and B in a magnetic field 
B. The change of the Lorentz force (a consequence of the changing 
sign of the carrier charge) creates a snake-shaped trajectory be- 
tween the two contacts, (b) False-color image of a device similar to 
the one studied, (c) The signature of the additional conduction is a 
reduction of resistance parallel to the junction: for four-terminal re- 
sistance measured parallel to the PNJ, Rodd (defined in the text) is 
negative (black trace), whereas the four-terminal resistance mea- 
sured perpendicular to the junction yields a postive Rodd (blue 
trace). 
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ers across the PNJ produces a change in the direction of 
the Lorentz force, causing classical trajectories to curve 
back toward the interface from both sides [Fig. 1(a)], sim- 
ilar to so-called "snake states" that propagate in 2D con- 
ductors along contour lines of an inhomogeneous mag- 
netic field [M]. Snake states in graphene have been con- 
sidered theoretically for inhomogeneous fields [151 US]? 
uniform fields at p-n interfaces [17 , and warped [18] and 
folded [19] graphene sheets. For the PNJ case [17], it 
was shown that classical trajectories similar to that of 
Fig. 1(a) can exist at the interface of p and n regions and 
that snake states continue to contribute to transport as B 
approaches zero. Related guiding of plasmons along PNJ 
interfaces has also been considered theoretically [20l |2T] . 

The main experimental observation of this Letter is il- 
lustrated in Fig. 1: while the presence of a PNJ increases 
resistance across a gate-controlled density interface, re- 
sistance along the same interface decreases when a PNJ 
is formed. Unlike the case of uniformly gated samples, re- 
sistances of gated PNJs are typically not symmetric with 
respect to the charge-neutrality point (CNP), since the 
top-gate voltage, Vt, controls density only on one side of 
the junction. Following Ref. [22j, we extract the contri- 
bution to resistance from just the PNJ by separating out 
the part of the four terminal resistance, R, that is odd 
relative to the CNP, Rodd(^T) = R{yT)-R(2V™^ -V^), 
where Kp^^ is the gate voltage that yields charge neu- 
trality under the top gate. When Vt yields a unipolar 
sample (with no PNJ), we take Rodd(^T) = 0. The well- 
known increase in resistance across a PNJ [6l [7] is ob- 
served by sourcing current at contact C with contact 1 
grounded, and measuring the voltage between contacts 
D and 2 [see Fig. 1(b) for geometry]. In Fig. 1(c), a PNJ 
is present for values of Vt to the right of the red dashed 
line, where Rodd is positive [blue trace in Fig. 1(c)], as 
expected in this configuration. In contrast, Rodd is neg- 
ative [black trace in Fig. 1(c)] for R measured along the 
PNJ, with current source at C, ground at D, and voltage 
measured between A and B. We interpret the negative 
value of Rodd in this configuration as the signature of an 
additional conduction channel at the p-n interface. 

Graphene fiakes made from mechanically exfoliated 
HOPG were placed on a degenerately doped Si substrate 
with a 300nm Si02 layer. Potential single layers were 
identified optically, and electrical contacts (5nm Ti/40nm 
Au) were subsequently patterned using electron-beam 
lithography and thermal evaporation. An unpatterned 
30nm AIO2 dielectric layer, produced by atomic layer 
deposition, with an NO2 functionalization layer [7 , was 
deposited on top of the graphene and the contacts. A top 
gate was patterned as shown in Fig. 1 using an aligned 
electron beam lithography step. The doped substrate is 
used as a back gate, controlling density globally; density 
under the top gate is set by voltage Vt (Region 1). Con- 
tacts straddled by the top gate [labeled A and B in Fig. 
1(b)] allow transport measurements along the PNJ. Mea- 



surements were made at temperature a of 4K using a 10 
nA current-bias measured using a lock-in amplifier. Two 
similar devices were measured, yielding similar results; 
data from one device is presented here. Magnetic fields 
up to 8T were applied perpendicular to the graphene 
plane. The device reported here shows QH signatures 
of single-layer graphene, i.e., conductance quantization 
at 2, 6, 10.... e^/h and has a CNP at Vb = 40V. Ah 
presented data were taken in the range -40V < Vb < 
20V (always p-type in Region 2, outside the top-gate re- 
gion) . Similar results were obtained (data not shown) for 
a smaller range of Vb > 40V. 

In Fig. 2, longitudinal resistance Sxx=Rcd,ab and 
transverse resistance Sxy=RAD,BC along the PNJ are 
compared to Rxx=Rcd,ef and Rxy=RcF,DE^ with all 
contacts beneath the top gate. {Rij,ki denotes a four 
terminal resistance with current applied at i and j and 
voltage measured at k and /.). Gate- voltage depen- 
dence of Rxx shows behavior typical of uniformly gated 
samples [23|: A fairly symmetric peak in Rxx around 
Vt^^^ - 3.6V identifies the CNP under the gate. For 
Vt< Vt^^^, the device is unipolar (p-p^); for Vt> 
V^rpCNP^ the device is bipolar, with a PNJ along an inter- 
face between contacts A and B. A 2D plot of Rxx{Vt^ B) 
[Fig. 2(a), inset] shows that while the resistance in- 
creases with increasing Vr^^^ does not depend on B. 
Sxx{Vt)^ measured along the PNJ, differs qualitatively 
from Rxx as well as from previously measured longitudi- 
nal resistances measured across PNJs Notably, in 
the bipolar regime, ^'^^^^(Vt) decreases by ^0.3kr^, pro- 
ducing a trace that is lower on the p-n side of the CNP, 
including at 5 = [Fig. 2(b), inset]. 

The Hall resistance in the top-gated region, Rxy{VT), is 
similar to data from single-gate graphene [23 [Fig 2(c)]: 
As the CNP is crossed, Rxy changes sign, indicating a 
change in carrier type {p ^ n) as a function of Vt- These 
curves are antisymmetric with respect to the CNP and 
B. On the other hand, the Hall resistance that involves 
the PNJ, Sxy^ is not antisymmetric with respect to the 
CNP; it is larger on the p-n side [Fig. 2(d)]. The dif- 
ference in Hall resistances, Sxy—Rxy is positive in the 
p-n regime throughout the measured range of 5, includ- 
ing zero [Fig. 2(e)]. A measurement similar to Sxy was 
previously used to study electron focusing in a 2D elec- 
tron gas [24]. Here, like in Ref. [2j, an increase in Sxy 
is indicative of enhanced transport between contacts A 
and B, as is observed here. It is important to note that 
the simultaneous increase in Sxy and decrease in Sxx are 
consistent with enhanced transport along the PNJ and 
rule out simpler explanations of the phenomena, like a 
position-dependent CNP within the sheet of graphene. 

In the quantum Hall regime {B > 4T), snake states 
also alter the nature of transport along a PNJ. Sxy in 
the unipolar regime, with the filling factor in Region 2 = 
-10, shows typical quantum Hall behavior for graphene, 
with plateaus at 1/2, 1/6 and 1/10 in units of h/e^ ^ and 
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FIG. 2: Measurements of the four quantities Rxx, Sxx, Rxy and Sxy as a function of Vt (black traces, offset intentional) for B between 
zb2T in 0.5T steps at a back-gate voltage of -20V. The measurement scheme is shown in the upper right inset (black contacts are the 
current leads and grey are the voltage leads) in each panel, (a) Rxx{Vt). The red dashed line locates the resistance maximum at all 
B fields and indicates the CNP for Region 1. Lower inset: RxxiVr, B) demonstrating that the CNP does not change over the B range 
explored, (b) SxxiVr)- A drop in resistance of 0.3 - 0.5kr2 is observed at the transition from p-p regime to p-n regime (red dashed line), 
indicating that an additional conduction channel has been introduced at the p-n interface. This resistance drop occurs for the entire B 
field range (lower inset), (c) RxyiVr). Hall resistance measurement where all the contacts are under the top gate reveals curves that are 
antisymmetric with respect to the CNP. In contrast, Sxy (d) has larger resistance on the p-n side of the CNP than Rxy, consistent with 
an additional conduction channel present at the p-n interface. This additional amount is quantified in (e) where a plot of Sxy— Rxy shows 
its largest value on the bipolar regime for all B fields (lower inset). 



depends only on the filling factor of Region 1, [black 
curve in Fig. 3(e), to the left of the red dashed line]. In 
the presence of a PNJ, the value of Sxy deviates from 
the typical values for graphene. When Region 1 is on a 
Hall plateau, i.e. for configurations shown schematically 
in Figs. 3(a,c), edge states equilibrate at contact A and 
are all ejected at the same voltage towards contact B. 
The values for Sxy are a result of edge state propagation 
and can be calculated via the Landauer-Buttiker formula, 
resulting in a close match with the experimental values 
of 2.3kQ and l.Gk^l for iyRi=2 and 6, respectively (see 
Supplementary Info). In addition to the modified Sxy re- 
sistances, a peak in resistance appears at Vt=5.5V at 8T 
[Fig. 3(e), indicated by black arrow] of magnitude ^2kQ. 
A smaller, but evident, peak is also visible at Vt=6.5 
V [smaller black arrow. Fig. 3(e)]. A plot SxyiVr^ B) 
shows that the stronger peak (large black arrow) moves 
linearly away from the CNP of Region 1 as 5 is increased 
[Fig. 3(d)] and follows the transition ^^^^1=2^6, suggest- 
ing that the peak is due to contributions to the resistance 
from pxx- This case is unlikely as the contribution from 
pxx in the p-p regime at the transition ujii=-6^-2 is 
not nearly as prevalent as this peak observed in the p-n 
regime. In between Landau levels, transport in Region 
1 is allowed to occur in the bulk, as is the case near 
B=0. In the same way that Sxy is enhanced at low fields 
in Fig. 2(d), it is enhanced here at the transition from 
^Ri=2^6 by the presence of a snake state at the interface 



[shown schematically in Fig. 3(b)]. Fig. 3(e) shows one- 
dimensional cuts of Fig. 3(d) starting at B=0 to B=8T 
in 2T steps, with each B- value curve shifted downward 
from the B=8T value for clarity. Here it is evident that 
the enhancement of Sxy at B=0 gradually evolves into 
the peak in the QH regime, suggesting that the low- and 
high-field phenomena have similar origins. Further, con- 
finement at the interface of the PNJ should be enhanced 
by a magnetic field, resulting in a stronger peak in Sxy, 
as is also evident in Fig. 3(e). 

Increasing Vb produces a voltage drop across the PNJ 
[AV=V(Region l)-V(Region 2), where both voltages are 
measured at the position of the peak in Sxy— Rxy] that 
increases, allowing for the transverse-electric-field depen- 
dence of the snake state to be investigated. The AV de- 
pendence of the peak resistance in Sxy— Rxy [from Fig. 
2(e)] is shown in Fig. 4, plotted for B between and 
2T, in 0.5T increments in the AV range of -23V to 83V 
(black lines are guides to the eye). It is found that the 
position of the peak moves linearly in the {Vt^Vb) space 
(data not shown) but decreases as the magnitude of AV 
is increased, i.e. as the electric field perpendicular to the 
junction increases. The change in resistance gets stronger 
for increases in the B field, changing ^0.2kQ at B=OT 
to —Ikl] at B=2T over the AV range shown here. 

An additional channel of conductance along the PNJ 
would alter the picture of charge transport at the charge- 
neutrality point, changing the value of (Jmin, the con- 
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FIG. 4: AV dependence of the difference Sx 
and 2T in 0.5T increments. Sx 
larger for all B fields. Solid black lines are guides to the eye. The 
decrease in resistance in this AV range increases as the perpendic- 
ular field is increased, rising from ~0.2kr2 at B=OT to ~lkr2 at 
B=2T. 
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FIG. 3: (a-c) Schematic of transport in the quantum Hall regime 
for a constant i>'r2=-^0, while Ufn progress from 2 (a), transition 
from 2 to 6 (b), and 6 (c). At the transition between filling factors 
2 and 6 (b), snake states enhance Sxy in a manner similar to Fig. 
2, with an amplitude that increases as B increases, (d) SxyiVr, B) 
reveals expected plateaus of 1/10, 1/6 and 1/2 /i/e^in the unipolar 
regime (Vt < 4V) and deviations from these values in the bipolar 
regime. An enhanced value of Sxy is observed at the transition 
between = 2 and 6 and is attributed to the presence of snake 
states at the p-n interface, (e) Constant-B-field cuts of (d) for 
B between and 8T in 2T steps. Cuts from to 6T are shifted 
downward from the B=8T cut for clarity. 



ductance at this point. Charge transport in disor- 
dered graphene samples has been studied experimen- 
tahy [H [27] and theoretical predictions have been made 
for (Jmin 1211 however consensus has yet to be reached. 
Taking into account the resistance of the PNJs and ap- 
propriate values for the size of density fluctuations, a 
value crmin"^2.5e^//i was obtained, 2 to 6 times lower 
than the experimentally reported values [3, 27 . The ad- 
ditional conduction along the interface could be a source 
conductance that brings this theoretical value closer to 
the experimental values. 
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